0-meson production in In-In collisions at £i a b=158A GeV: evidence for relics of a 

thermal phase. 



OS 

o 
o 

(N 

<D 

in 

(N 



E. Santini, H. Petersen, M. Bleicher 
Institut fur Theoretische Physik, Goethe-Universitat, 
Max-von-Laue-Str. 1, D-60438 Frankfurt am Main, Germany 

(Dated: September 25, 2009) 

Yields and transverse mass distributions of the ^-mesons reconstructed in the <f> ~^ 
channel in In+In collisions at E\ a b=158A GeV are calculated within an integrated Boltz- 
mann+hydrodynamics hybrid approach based on the Ultrarelativistic Quantum Molecular Dynamics 
(UrQMD) transport model with an intermediate hydrodynamic stage. The analysis is performed 
for various centralities and a comparison with the corresponding NA60 data in the muon channel is 
presented. We find that the hybrid model, that embeds an intermediate locally equilibrated phase 
subsequently mapped into the transport dynamics according to thermal phase-space distributions, 
gives a good description of the experimental data, both in yield and slope. On the contrary, the 
pure transport model calculations tend to fail in catching the general properties of the <j> meson 
production: not only the yield, but also the slope of the bit spectra, very poorly compare with the 
experimental observations. 
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I. INTRODUCTION 

The production of <fi mesons is considered to be one 
of the key observables to probe the state of matter pro- 
duced in relativistic heavy-ion collisions. Strangeness en- 
hancement in relativistic nucleus- nucleus collisions com- 
pared to nucleon-nucleon collisions has been originally 
suggested as a possible signal for the formation of a de- 
confined plasma of quarks and gluons during the initial 
state of the reaction [l], 0, Q • The dominant production 
of ss pairs via gluon-gluon interaction in the plasma may 
result in an enhanced number of strange and multistrangc 
particles produced after hadronization; in particular, free 
ss pairs would coalesce to form <j> mesons [jj, whereas 
their production in pp collisions is suppressed according 
to the Okubo-Zweig-Iizuka rule [4|, la, |6|. It is moreover 
expected that <j> mesons decouple from the rest of the 
system earlier than other non-strange hadrons. At RHIC 
energies, an early decoupling of the (f> from the hadronic 
rescattering dynamics was found in Refs. 0] and Q. 
Similar conclusions were obtained with the RQMD cas- 
cade in Q for the f2 baryons at SPS energies. 

The dominant hadronic decay of the <j> meson is <fi ~~ * 
KK; additionally, being a neutral vector meson, the <j> 
contributes to dilepton production via the direct decays 
<t> — * e + e~ and — > fi + fj," . Phi meson production was in- 
vestigated extensively at the SPS by several experiments 
in the kaon (NA49_Ip [H Ql and CERES (3), di- 
electron (CERES [lj), and dimuon (NA50 [3 and NA60 
(l6| ) channels. The reconstruction ofin-matter (j) - * KK 
decays, however, might be partially prevented by kaon 
absorption and rescattering [l?], HI, HH and a priori a 
careful investigation of kaon final state interactions can- 
not be avoided in a quantitative comparison to experi- 
mental data. On the contrary, the dilepton channel is 



not affected by such a shortcoming: dileptons leave their 
production site essentially undistorted, and a compari- 
son of model calculations to measurements in the dilep- 
ton channel is surely more straightforward. Recently, 
high statistics measurements of <f> — > fi + fj,~ meson pro- 
duction in In-In collisions at £'i a b=158A GeV have been 
presented by the NA60 Collaboration 16]. This reaction 
will be addressed in the present work. 

To investigate <f> — » ^t + /i~ production in In-In colli- 
sions at £ ; i a b=158A GeV we employ an integrated Boltz- 
mann+hydrodynamics hybrid approach based on the Ul- 
trarelativistic Quantum Molecular Dynamics (UrQMD) 
transport model with an intermediate hydrodynamic 
stage. In this approach, initial conditions and contin- 
uous decoupling up to freeze-out are treated within a full 
non-equilibrium transport approach, whereas hydrody- 
namics is used to describe the intermediate equilibrated 
phase. This allows to reduce the parameters for the ini- 
tial conditions and provides a consistent freeze-out de- 
scription. Moreover, by comparing the hybrid approach 
to pure transport calculations, we are able to directly 
investigate the consequences that a dynamical approach 
involving local thermal and chemical equilibrium and one 
based on full non-equilibrium dynamics have on <p meson 
production. 

The paper is structured in the following way: In Sec. 
HIl we briefly discuss the hybrid model and present the 
procedure used to evaluate the <fi — > /i + [i~ emission and, 
consequently, reconstruct the <f> meson. In Sec. IIIII we 
present calculations for </>-meson transverse mass spectra 
as a function of centrality. The various contributions to 
the spectra associated with different stages of the reac- 
tion dynamics are shown and a comparison to the NA60 
data is presented. Finally, a summary and conclusions 
are given in Sec. IIV1 
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II. THE MODEL 



A. The hybrid approach 



To simulate the dynamics of the In+In collisions we 
employ a transport approach with an embedded three- 
dimensional ideal relativistic one fluid evolution for the 
hot and dense stage of the reaction based on the UrQMD 
model [20l | . The present hybrid approach has been exten- 
sively described in Ref. [20j. Here, we limit ourselves to 
briefl y d escribe its main features and refer the reader to 
Ref. [20] for details. 

UrQMD HH, S3, [H is a hadronic transport ap- 
proach which simulates multiple interactions of ingoing 
and newly produced particles, the excitation and frag- 
mentation of color strings and the formation and de- 
cay of hadronic resonances. The coupling between the 
UrQMD initial state and the hydrodynamical evolution 
proceeds when the two Lorentz-contracted nuclei have 
passed through each other. Here, the spectators con- 
tinue to propagate in the cascade and all other hadrons 
are mapped to the hydrodynamic grid. This treatment is 
especially important for non-central collisions which are 
also studied in the present work. Event-by-event fluc- 
tuations are directly taken into account via initial con- 
ditions generated by the primary collisions and string 
fragmentations in the microscopic UrQMD model. This 
leads to non-trivial velocity and energy density distribu- 
tions for the hydrodynamical initial conditions [13, H3| • 
Subsequently, a full (3+1) dimensional ideal hydrody- 
namic evolution is performed using the SHASTA algo- 
rithm [26l . |27| . The hydrodynamic evolution is gradu- 
ally merged into the hadronic cascade: to mimic an iso- 
eigentime hypersurface, full transverse slices, of thickness 
Az = 0.2fm, are transformed to particles whenever in 
all cells of each individual slice the energy density drops 
below five times the ground state energy density. The 
employment of such gradual transition allows to obtain a 
rapidity independent transition temperature without ar- 
tificial time dilatation effects [28[ and has been explored 
in detail in various recent works [28|, H(J HH devoted 
to SPS conditions. When merging, the hydrodynamic 
fields are transformed to particle degrees of freedom via 
the Cooper-Frye equation in the computational frame. 
The created particles proceed in their evolution in the 
hadronic cascade where rescatterings and final decays oc- 
cur until all interactions cease and the system decouples. 

An input for the hydrodynamical calculation is the 
equation of state (EoS). In this work we employ a hadron 
gas equation of state, describing a non-interacting gas of 
free hadrons [32|. Included here are all reliably known 
hadrons with masses up to sa 2 GeV, which is equivalent 
to the active degrees of freedom of the UrQMD model. 



B. Reconstruction of the phi meson in the dimuon 
channel 

The reconstruction of the <j> meson from the dimuon 
channel requires the evaluation of the <j) — > emis- 
sion. This is calculated perturbatively in the evolution 
stage that precedes or follows the hydrodynamical phase, 
and from thermal rates in the latter. In the following, the 
terms pre-equilibrium/pre- hydro (post-equilibrium/post- 
hydro), will be used to indicate the stage preced- 
ing (following) the mapping UrQMD^hydrodynamical 
(hydrodynamical— >UrQMD) evolution description. Note 
that in the pre- and post-hydro stages the particles are 
the explicit degree of freedom and their interactions are 
explicitly treated within the cascade transport approach. 

1. Reconstruction from pre- equilibrium and 
post-equilibrium emission 

Given the number of /i + fi pairs emitted in 

<f> — ► decays, the number iVT ec of <f> meson recon- 

structed in the dimuon channel is 

where BR^ + ^_ is the cf> — > branching ratio. The 

latter is given by the ratio between the dimuon and the 
total width of the tf> meson, i.e. BR^ + ^_ = r^ + ^_/rf ot . 
Thus, 

N reo = N <k 4r^- (2) 
H+ft- 

In the pre-equilibrium and post-equilibrium phase 
dimuon emission from the </>-meson can be calculated per- 
turbatively as 




where now indicates the number of (f> mesons present, 
in some stage of the evolution, in the system and Ti(n) 
{jf{n)) are the times at which the n-th </> meson appeared 
in (disappeared from) the system and are evaluated in the 
meson rest-frame. This perturbative method is known 
as time integration method or "shining method" and has 
long been applied in the transport description of dilepton 
emission (see e.g. [H, (Mill)- Combining Eqs.Q-© 
one has: 




which expresses the fact that the number of reconstructed 
mesons is proportional to the typical life-time of the me- 
son in the system. 
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Some considerations are now in order. During the pre- 
cquilibrium phase, Ti(n) coincides with the time at which 
a 4> meson is typically produced from a nucleon-nucleon 
scattering. With exception of some very rare almost in- 
stantaneous interaction, 77(71) typically coincides with 
the time at which the hydrodynamical phase starts. In 
other words, the probability of dimuon emission in the 
pre-hydro stage is evaluated up to the moment 4> mesons 
are merged into the hydrodynamical phase. From then 
on, (j> emission will be treated as thermal. 

In the post-equilibrium phase, Tj(n.) coincides either 
with the time at which the transition from the hydro- 
dynamical to the transport description is performed (for 
those (j) mesons produced via the Cooper-Frye equation) 
or with the time at which a <f> meson is produced from 
hadronic interactions still occurring in the cascade phase, 
as result of the fact that the whole system is not yet com- 
pletely decoupled; 77(71) is, in this case, the time at which 
the 4> meson decays or eventually rescatters. 



2. Reconstruction from equilibrated thermal emission 

Thermal dimuon emission from the <j> meson can 
be expected to be significantly smaller than the post- 
equilibrium emission, due to the fact that the lifetime of 
the fireball (7-10 fm) is much smaller than the <j> (vac- 
uum) lifetime of ^44 fm. To determine the thermal 
dimuon production rate from 0-meson decays we observe 
that the mass of the <f> meson (to0=1.O19 GeV) is larger 
than the typical local temperature of the thermalized fire- 
ball, thus the particle number distribution function can 
be reasonably evaluated in Boltzmann approximation. 
Moreover, the <j> meson being a very narrow resonance 
(rfot = 0.00426 GeV and rf ot /m w 0.4% ), for simplic- 
ity we neglect its small width and approximate the mass 
distribution of the meson with a 8(m 2 — m 2 ,) function 
(pole approximation). 

In the Boltzmann and pole approximation the <fr par- 
ticle phase space distribution function is given by 



dN A 3 
{q ;T[x,t)) = 



d 3 x d 3 q 



(27T) 



,-<R,/T(x,t) 



(•5) 



where T(x, t) is the local temperature, (go> q) the meson 
4-momentum in the thermal frame and the dependence of 
the temperature from the (discretized) space-time point 
of the (3+1) grid has been explicitly indicated. The num- 
ber of dimuons produced per unit phase space volume 
and unit time from <j> — > fi + fj,~ decays is 



d 3 x d 3 q 70 (q ) 



(6) 



Here T q 



is the dimuon width as defined in the meson 



rest-frame and 7^ = qo/m^ is the Lorentz factor trans- 
forming from the meson to the thermal rest frame. The 
number of dimuons emitted per unit space-time volume 



can be obtained by integrating Eq. (J6]) over momentum. 
The integration can be performed analytically and one 
finds: 



^7-(T(x,0) 



3 

2^2 



miT(x,t) K\ 



T(x,t) 



(7) 

where K\ is the modified Bessel function of first order. 
The contribution of a single cell to the dimuon emission 
in the time step At is therefore: 



'cell 



dN M+ll - 
d A x 



\T ce \\) Vcell At. 



(8) 



The cell contribution to the dimuon emission is calculated 
according to the above expression. To compare with ex- 
perimental data, dimuon momenta in the cm. frame are 
then generated with a Monte Carlo procedure according 
to the distribution function 



d 3 x d 3 p po 



JB 



T r <,n 



(9) 



with u v the fluid cell 4-velocity and fs the boson distri- 
bution function. 

The total dimuon yield from thermal (f> mesons is ob- 
tained summing the above expression over all fluid cells 
and all time steps of the (3+1) grid that are spanned 
by the system during the hydrodynamical evolution un- 
til the transition criterium is reached. Let us express this 
summation symbolically as: 



= f d 4 / 3 2 T( t) K ( _J^±_] r <t> 

'10; 

The corresponding number of <fi mesons reconstructed in 
the dimuon channel is then: 



JJL mlT(x,t) Kx ( 
\2tt 2 * V ' ' \T(x,t) 



NjT = / d^x 

J VCR 4 

(11) 

Before proceeding, we would like to make two consider- 
ations. The first one concerns an implication of the EoS 
used for the hydrodynamical evolution. As previously 
stated, a hadronic EoS has been used in the present work. 
This implies that dimuon emission from the <p meson is 
here evaluated during the whole hydrodynamical phase. 
In general, however, if a phase transition from quark- 
gluon to hadronic matter occurs, the hadronic thermal 
rate would be only a fraction of the total thermal rate. 
In this sense, the results presented in the next section 
on thermal dimuon emission from the 4> meson should 
be regarded as an upper limit of the possible amount of 
dimuons indeed produced by (f> mesons during the high 
temperature/high density stage of the heavy-ion colli- 
sion. 

The second consideration regards the validity of the 
pole approximation used. This approximation is justi- 
fied as long as the <f> meson maintains its properties of 
narrow resonance. In medium, the <p meson is expected 
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to broaden, as suggested by hadronic many body cal- 
culations. If the amount of broadening is significant 
or/and the spectral function of the meson develops a 
complex "structured" shape in the medium, the pole ap- 
proximation may loose validity. On a quantitative base, 
there is no general consensus on the specific amount of 
broadening of the <j) meson spectral function to be ex- 
pected in a high temperature/high density environment. 
Early evaluations of medium effects based on hadronic 
rescattering at finite temperature have indicated quite 
moderate ch ang es of both the <f> meson mass and width 
[3(1 l37l . [H, |39| ] . In a subsequent calculation, collision 
rates in a meson gas have been estimated to amount to a 
broadening by ~ 20 MeV at T=150 MeV [ZtJ. The dress- 
ing of the kaon cloud is presumably the main effect for <j> 
modifications in nuclear matter, increa sing its width by 
~ 25 MeV at normal nuclear density [4JJ]. In hot and 
baryon poor hadronic matter the in-medium properties 
of the <p have been schematically explored in Ref. [42| : 
the meson was found to retain its resonance structure and 
an in-medium width of ~ 32 MeV at (7Vuat)=(180,27) 
MeV was estimated. In a recent work [4j], the spectral 
density of the cj> meson in a hot bath of nucleons and pi- 
ons has been microscopically calculated from the forward 
scattering amplitude in a two component approach. The 
authors found a considerable broadening of the meson 
width, e.g., r^ lcd ~ 100 MeV at saturation density and 
temperature T = 150 MeV. 



In fact, there is no visible evidence for a strong in- 
medium scenario for the <p meson from the NA60 data. 
For all the analyzed centrality bins, the measured in- 
variant mass distribution can be described in terms of 
a vacuum-like spectral function and the extracted val- 
ues for the mass and the width are compatible with the 
PDG values and independent of centrality. Of course, 
one should remind that the extraction of the in-medium 
modified component of the total dimuon emission from 
the experimental data is a non-trivial task, since this 
component would lie under the large unmodified peak 
produced from the decays occurring at the freeze-out. 
Certainly, the study of in-medium modifications of the <j> 
meson properties is in itself an interesting research topic. 
However, in the present work will not address this is- 
sue and assume that also in-medium the resonance main- 
tains its narrow width, so that the pole approximation 
is still valid. For the meson pole mass m^, moreover, 
the vacuum value will be used. In any case, as we will 
show below, from our analysis it emerges that the amount 
of cf)'s from thermal emission is by far smaller than the 
abundance from the cascade part, so that an eventual in- 
medium modification of the thermal rate will most likely 
not alter the results on the total transverse mass spectra 
discussed in the next section. 



III. RESULTS 

In this section we investigate the relative abundances 
of <j> — > production in the various stages of the 

system evolution and present results for 4> meson invari- 
ant transverse mass spectra as a function of the collision 
centrality. Calculations for cj> — > production in 

In-In collisions at E\ a \y—158A GeV have been performed 
for 5 centrality classes. In agreement with the treatment 
of the experimental data, each class was identified by 
the range of charged particle multiplicity dN c h /drj in the 
pseudorapidity interval 3< rj <4. The relation between 
a specific range of dN c h/dr) and the corresponding cen- 
trality bin was specified by the NA60 collaboration and 
can be found in Table 1 of Ref. [l(| . The correspondence 
between the 5 ranges of dN c h/drj and 5 ranges of impact 
parameters was obtained by the analysis of the charged 
particles obtained within the UrQMD hybrid model as a 
function of the impact parameter selected in the Monte 
Carlo simulation. 



A. Relative abundances in the various evolution 
stages and hit spectra 

First, let us investigate the relation between the 
amount of <p mesons reconstructed from cj> — > fj, + fi~ and 
the stage of the evolution probed by the dimuon emission. 
We start by discussing separately and in detail results for 
the most central bin. Later, an analogous analysis is pre- 
sented as function of centrality. 

The contribution to the transverse mass spectra of 
those 4> mesons emitting during the hydrodynamical and 
the cascade stage are separately shown in FigfTJ The 
emission from the hydrodynamic stage is found almost 
one order of magnitude smaller than the emission from 
the cascade. As already mentioned, this is due to the 
relative smallness of the duration of the hydrodynam- 
ical phase when compared to the cascade phase. Thus, 
the dominant contribution comes from the cascade stage. 
This stage contains pre- and post-equilibrium emission. 
However, as shown in Figf3J the pre-equilibrium emis- 
sion is negligible (two orders of magnitude smaller than 
the post-equilibrium emission) and will not be discussed 
further. The post-equilibrium emission can be divided 
in two categories: (i) the emission from (f> particles pro- 
duced via Cooper-Frye at the transition point and (ii) 
the emission from <fi particles produced during the cas- 
cade (Fig. [2]). In the first case, the particles have a 
momentum distribution that reflects the thermal prop- 
erties of the transition point, although their dimuon de- 
cay occurs later in a non-thermal environment. In this 
sense, this copious emission, though not specifically ther- 
mal (i.e. not described by thermal rate equations) still 
carries information about the preceding thermal phase. 
It represents the most direct "remain" of the thermally 
equilibrated phase previously experienced by the system. 
The second contribution, on the contrary, can be labeled 
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FIG. 1: Transverse mass distributions of the <f> meson in cen- 
tral indium-indium collisions. The <f> is reconstructed in the 
dimuon channel (see text). Double-dotted-dashed line: con- 
tribution to the (j> — > /j> + n~ production of the hydrodynamical 
stage. Dashed line: contribution of the cascade stage. Full 
line: total <j) — > fi + n~ production. 

as a "purely cascade" one. This is the contribution of 
(j> particles produced in the non-equilibrium environment 
on the way to final decoupling. This second contribution 
is characterized by steeper transverse mass spectra, the 
shape of the total spectra is found to be composed by the 
interplay of both emissions. 

It is instructive to compare the hybrid model calcula- 
tions to the pure cascade calculations, in which no as- 
sumption of an intermediate equilibrium phase is made. 
The comparison is presented in Fig. [3] (top) together 
with the experimental data. As one can see, the absence 
of an intermediate thermal phase results in a steepening 
of the transverse mass spectra not supported by the data. 
Moreover, the pure cascade calculation strongly underes- 
timates the <p meson yield, a feature already emerged in 
recent independent investigations [l3| performed in rela- 
tion to measurements obtained by the NA49 experiment. 
There, it was found that a statistical hadron gas model 
with undersaturation of strangeness [44j could account 
for the measured yields [13| . 

From this first analysis, we can conclude that, despite 
the smallness of the specifically thermal <f> — > /J + /i~ emis- 
sion, the presence of the thermal phase is essential in or- 
der to obtain an appropriate yield and slope of the Tot 
distributions from the cascade emission. In this sense, 
we speak about the presence of "thermal relics" . 

An analogous analysis has been performed for the fur- 
ther 4 centrality bins and results are presented in Fig[4] 
and FigO For all centrality classes the thermal rate from 
the hydrodynamic evolution is found to be much smaller 
than the cascade emission. With decreasing centrality, 



FIG. 2: Decomposition of the post-equilibrium (j> — > 
production in: (i) emission from (f> mesons which are merged 
into the cascade at the transition point via the Cooper- 
Frye equation and emit in the cascade stage of the evolution 
(dashed line) ; (ii) emission from <f> mesons which are produced 
and emit in the cascade stage (dotted-dashed line). The pre- 
equilibrium emission is denoted by the dotted line. The full 
line represents the total cascade emission. Due to the small- 
ness of the pre-equilibrium emission the latter practically co- 
incides with the total post-equilibrium emission. 
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FIG. 3: Transverse mass distributions of the (j> meson in cen- 
tral (top) and peripheral (bottom) indium-indium collisions. 
The hybrid model calculation (full line) is compared to the 
pure UrQMD transport calculation (dashed line) and to ex- 
perimental data [16| . Bin- widths which coincide with the ones 
of the experimental data have been here used. 
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FIG. 4: Same as in Fig. [T] but for different centrality classes. 
Double-dotted-dashed line: contribution to the <j> — * fi + yT 
production of the hydrodynamical stage. Dashed line: con- 
tribution of the cascade stage. 



10- 



CS- 10 1 

'> 

a io- 1 

5 -3 

£ 10 3 

? 

B 10 1 



10" 



10" 



r In-ln n 
r 95<dN ch /dn<16(h 

I , I , I , I , I , I , I , I , I , = 


r In-ln -j 

r 50<dN ch /dr|<95 -j 

-,1,1,1,1,1, ,\, : 


r In-ln i 

r 25<dN ch /dT]<50 -j 

I , I , I , I , I , I , I , l \ r 


r In-ln i 
r 4<dN ch /dr|<25 -j 




10 



In-ln 

" (|)-meson reconstructed from § -> \ 

, i , , , , i , , , , i , , , , i 



0.5 1 1.5 2 2.5 
rrym [GeV] 



0.4 0.8 1.2 1.6 0.4 0.8 1.2 1.6 2 
nynyGeV] rryrryGeV] 



FIG. 6: Transverse mass distributions of the (f> meson in 
indium-indium collisions as a function of centrality; from top 
to bottom: central to peripheral spectra. The hybrid model 
calculations are compared with NA60 data [l6j . Bin-widths 
which coincide with the ones of the experimental data have 
been here used. 



latccl within the hybrid model for the 5 centrality classes 
considered are compared to NA60 data in Fig. [SI The 
hybrid model can account pretty well for both slope and 
yield in the first 4 centrality classes. Small deviations 
are observed for the utmost peripheral bin where, most 
likely, this kind of hybrid models are already at the limit 
of applicability. They rely on the assumption that an 
equilibrium phase is indeed reached, which for very pe- 
ripheral reactions is at least questionable [45[ . 

Pure transport calculations fail however too in describ- 
ing this very peripheral reaction (see Fig. [31 bottom), 
though the comparison with experimental data suggests 
that deviations of the resulting slope from the measured 
one are much smaller than the in the case of central col- 
lisions. 



FIG. 5: Same as in Fig. [2] but for different centrality classes. 
Dashed line: emission from (j> mesons which are merged into 
the cascade at the transition point via the Cooper- Frye equa- 
tion and emit in the cascade stage of the evolution; Dotted- 
dashed line: emission from (j> mesons which are produced and 
emit in the cascade stage. 

the pure cascade emission (dash-dotted line) becomes less 
and less important and the spectra is more and more de- 
termined by the emission from those mesons emerging 
from the thermal stage into the cascade at the transition 
point. 

Finally, the (j) meson transverse mass spectra calcu- 



IV. SUMMARY AND CONCLUSIONS 

In this work we employed an integrated Boltz- 
mann+hydrodynamics hybrid approach based on the 
Ultrarelativistic Quantum Molecular Dynamics trans- 
port model with an intermediate hydrodynamic stage 
to analyze ^-meson production in In+In collisions at 
-Eiab=158A GeV from its reconstruction in the <fi — > 
channel. We find that the hybrid model fairly de- 
scribes </> yields and transversal mass spectra at various 
collision centralities. In particular, the analysis points 
out that the underlying assumption of the existence of an 
intermediate equilibrated phase seems to play an essen- 
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tial role in order to catch the main aspects of the physics 
emerging from the <f> dimuon emission at top SPS ener- 
gies. 
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